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Abstract 

A search for the supersymmetric partner of the goldstino, the sgoldstino S, at 
LEP2 is presented. The production S7 followed by S decay into two gluons or 
into two photons was studied at 189 - 202 GeV LEP centre-of-mass energies. No 
evidence for the S production was found and limits on the S mass corresponding 
to different theory parameters are given. 
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1 Introduction 

In the Supersymmetric extension of the Standard Model, once Supersymmetry is spon- 
taneously broken the gravitino G can acquire a mass by absorbing the degrees of free- 
dom of the goldstino. This mechanism is analogous to the spontaneous breaking of the 
electro-weak symmetry in the Standard Model, where the Z and W bosons acquire mass 
by absorbing the goldstone bosons. 

A light gravitino as predicted by some supersymmetric rnodels ^ has been searched 
for at LEP and Tevatron experiments Limits on the G mass allow lower limits on 

the supersymmetry-breaking scale a/F to be inferred. 

Recently it has been pointed out [Q that an appropriate theory must contain also 
the supersymmetric partner of the goldstino, called the sgoldstino, which could be mas- 
sive. The production of this particle may be relevant at present LEP energies if the 
supersymmetry-breaking scale and the sgoldstino mass are not too large. In the minimal 
R-parity-conserving model, as considered in the effective theory at the weak scale 
contains two neutral scalar states: the 5* which is CP-even, and the P which is CP-odd. 
As sgoldstinos have even R parity, they are not necessarily produced in pairs and their 
decay chains do not necessarily contain an LSP (Lightest Supersymmetric Particle). The 
phenomenology of these two particles is similar. The following formulae and results will 
be expressed for the S state but are valid also for the P particle. 

At LEP 2, one of the most interesting production channels is the process e'^e" — > 6*7, 
which depends on the S mass nis and on ^/F : 

'e-e- -^S,) = ^(l-^]\l + cos^e), (1) 



s 



where 6 is the scattering angle in the centre-of-mass and 



with Ve = sin^Ow — 1/4, ag = 1/4 and gz = e/ {sin9wcos9w) ■ and M^z are related 
to the diagonal mass terms for the f/(l)y and SU{2)l gauginos Mi and M2: 

M^^ = Micos^dw + M2sin'^9w, M^z = (M2 - Mi)sin9wcos9w (3) 
The most relevant S decay modes are S" — 77 and S ^ gg with 

ns ^ 77) ^ ^ (4) 



327rF2 



and 



r(5 - 99) = (5) 

where M3 is the gluino mass. The corresponding branching ratios depend on Mi, M2 
and M3, and the total width is F ~ r(S' — 77) + r(S' — gg). In this letter two sets for 
these parameters as suggested in [|] are considered and listed in Table |l]. 

For a large interval of the parameter space the total width is small (below a few 
GeV/c^), except for the region with small a/F where the production cross section is also 
expected to be very large. 

The two decay channels considered produce events with very different topologies: 
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Ml 


Ms 


M3 


B.R. ^ ^ 77 


B.R. S^gg 


1) 


200 


300 


400 


4% 


96% 


2) 


350 


350 


350 


11% 


89% 



Table 1: Two choices for the gaugino mass parameters (in GeV/c^) relevant for the 
sgoldstino production and decay and the corresponding Branching Ratios of the two 
considered channels. 

1. S* — i> 77 gives rise to events with three high energy photons, one of which is expected 

2 

to be monochromatic with energy = -^f- for the large fraction of the parameter 
space where S has a neghgible width. Despite its lower branching ratio (4 and 
11% for the two sets of Table |I], respectively), this final state is worth investigating 
because the main background source is the QED process e+e~ — > 77(7), which is 
expected to be small if photons in the forward region are discarded. 

2. S ^ gg gives rise to events with one monochromatic photon (except for the region 
with small \fF ) and two jets. An irreducible background from e^e~ qq'y events 
is associated to this topology. Therefore the signal must be searched for as an excess 
of events over the background expectation for every mass hypothesis. 

This letter describes the results obtained with the DELPHI detector at LEP centre- 
of-mass energies of 189, 192, 196, 200 and 202 GeV, corresponding to a total integrated 
luminosity of about 380 pb^^. 

2 Apparatus 

A detailed description of the DELPHI detector can be found in 0. The present 
analysis was mainly based on the measurement of the electromagnetic energy clusters 
in the barrel electromagnetic calorimeter, the High density Projection Chamber (HPC), 
and in the Forward ElectroMagnetic Calorimeter (FEMC), as well as on the capability 
of reconstructing charged particle tracks using the tracking devices: the Vertex Detector 
(VD), the Inner Detector (ID), the Time Projection Chamber (TPC), the Outer Detector 
(OD) and the forward chambers (FCA and FCB). The Vertex Detector extends its 
coverage down to 10.5° in polar angle 6. An electromagnetic calorimeter (STIC) was used 
to measure the luminosity. 

The barrel and the forward electromagnetic energy triggers were based on data from 
the HPC and the FEMC respectively. The calorimetric trigger efficiency for e+e~ — > 
77(7) was estimated with samples of Bhabha e~^e~ —>■ e"'"e~(7) events. This was done by 
counting how often the electromagnetic trigger was fired by an electron which had been 
triggered by an independent track trigger. In events with more than two photons, as well 
as in events with photons and charged particle tracks, the trigger efficiency was better 
than 99%. 

3 Event selection and analysis 

The 1998 data were taken at = 188.6 GeV, and the 1999 data at 191.6, 195.5, 
199.5 and 201.6 GeV. The integrated luminosities obtained requiring the HPC, FEMC, 
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TPC and VD to be operational were 155.1 pb~ , 25.1 pb~ , 76.2 pb^ , 83.1 pb^ and 
40.1 pb~^ respectively for the five centre-of-mass energies. 

Monte-Carlo generated events for the same centre-of-mass energies were processed 
through the full DELPHI simulation ^ and the same reconstruction chain as real data. 



3.1 5 ^ 77 channel 

Events were selected as 777 candidates if they had: 

• at least two electromagnetic energy clusters with 0.219 < E / y/s < 0.713; 

• at least one additional cluster with E > 5 GeV and no more than two additional 
clusters, of which the second one (if present) had E < 5 GeV; 

• the two most energetic electromagnetic clusters in the HFC region 42° < 6 < 89° or 
in the FEMC region 25° < ^ < 32.4°; 

• the third cluster in the region 42° < 9 01 20° < 9 < 35°; 

• no hits in two of the three Vertex Detector layers within ±2° in azimuthal angle (j) 
of the line from the mean beam crossing point to any electromagnetic cluster. 

Further, two hemispheres were defined by a plane orthogonal to the direction of the 
most energetic cluster. One hemisphere was required to have no charged particle detected 
in the barrel region of the detector with momentum above 1 GeV/c extrapolating to 
within 5 cm of the mean beam crossing point. The requirement was strengthened, to 
suppress the large e~^e~ background further, by demanding that both hemispheres have 
no such particle detected by the TFC with 9 < 35°. 

The events selected have a three-body final state kinematics if no significant additional 
radiation is lost in the detector (mainly initial state radiation lost along the beam pipe). 
A simple way to check if an event is, within a reasonable approximation, a three-body 
final state, is to look at the distribution of the quantity A = |5i2| + l^isl + 1^23 1 where 
5ij is the angle between the particle i and j (Fig. |l|). In a three-body final state, the 
particles lie in a plane and therefore A should be 360°. If only the events with A > 358° 
are accepted, the energies of the particles can be determined with very good precision 
from their measured directions: 

El = Vs — - — ; E2 = Vs — — ; ^3 = vs — — (6) 

with S = sinSi2+sin6is+sin623. The error on the energy evaluation was further minimised 
by requiring mm(5i2, 5i3, ^23) > 2°. 

In 5*7 events, the 5* decay products are isotropically distributed in the S centre-of- 
mass system. The distribution of cosa, where a is the angle between the S direction 
(opposite to the prompt photon) and the direction of one of the two S decay products, 
in the S centre-of-mass system, should therefore be flat. On the other hand, in the QED 
background |cosa| peaks at 1. Therefore, out of the three combinations present in each 
event, only those giving \cosa\ < 0.9 were accepted. 

The numbers of selected events, each giving up to three combinations, are listed with 
the expected background in Table 0. No significant background was found except for the 
QED process e~^e~ — > 77(7). 

The acceptance for an S-y signal produced according to (P after the described polar 
angle cuts was (51 ± 2)%. The dependence on ms from 10 to 190 GeV/c^ was contained 
within the error quoted. The selection efficiency inside the acceptance region was evalu- 
ated by means of the QED background events generated according to M, The efficiency 
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Figure 1: A — 180° for the 777 candidates (points) and the QED e~^e~ 77(7) simulated 
sample (histogram). The cuts on ™n(5i2, ^13, ^23) > 2° and on cosa are not applied in 
this figure. 

was independent, within the errors, of the photon polar angle. Its average value was 
(76.6 ±2.5)%. 

The energy resolution obtained from (^ was also evaluated using simulated QED 
events as shown in Fig.|^. It was better than 0.5 % in the whole photon energy range: a 
fit with two Gaussians gave two resolution components with (Xi = 0.12 GeV and a2 = 0.35 
GeV with about equal frequencies. 

The second component was introduced to describe the tails originating from photons 
detected near the calorimeter dead regions. 



3.2 S ^ gg channel 

This channel is expected to give rise to a final state with one photon and two jets. An 
event was selected as a 'ygg candidate if it had: 

• an electromagnetic energy cluster identified as a photon with E > 5 GeV and 9 > 
20°; 

• no electromagnetic cluster with ^ < 5°; 

• total multiplicity greater than 10; 

• charged particle multiplicity greater than 5; 

• \/J2i=iipl + Pl)i > 0.12 X yi, where n is the total multiplicity; 

• the sum of the absolute values of all particle momenta along the thrust axis greater 
than 0.20 x ^; 
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Figure 2: The energy resolution AE for the photons of the 777 candidates in the QED 
e^e~ — i> 77(7) simulated sample. The photon energy was obtained using d). A fit with 
two Gaussians gave two resolution components cxi = 0.12 GeV and a2 = 0.35 GeV with 
approximately equal frequencies. 



• either an electromagnetic cluster with E < 0.45 x ^/s, or a total multiplicity greater 
than 16 if the cluster energy is greater than 0.45 x ^/s ; 

• \cos{9p)\ < 0.995, where 9p is the polar angle of the missing momentum; 

• visible energy greater than 0.60 x ^/s^, 

• I cosa| < 0.9; 

• A greater than 350°. 

The events were reconstructed forcing all particles but the photon into a 2-jet topology 
using the DURHAM algorithm. Events were removed if ycut > 0.02 and if the angle 
between the photon and the nearest jet was less than 10°. If the event contained more 
than one photon candidate, the most energetic one was considered as the one produced in 
e'^e" S'-y. In addition, the jets were required to be incompatible with the bb hypothesis 
by requiring the combined btag of the events to be less than zero 

As in the 777 selection, the events obtained after this selection are three-body final 
state events in the absence of additional lost radiation. Therefore all the kinematic 
constraints described in the previous subsection were also applied here. In this case, 
however, as jet directions are less precisely determined than photons directions, the cut 
in A was less stringent and the resolution for the reconstructed photon energy was poorer: 
a two-Gaussian fit gave cxi = 1.2 GeV (55% of the area) and (72 = 4.1 GeV. 

The polar angle acceptance for an 5*7 signal produced according to (P was (76 ± 2%) 
and almost independent of m^. The selection efficiency inside the acceptance region 
was evaluated using the qq'-f background events generated with PYTHIA [O, processed 



through the full DELPHI analysis chain and re-weighted according to the background 
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and signal photon polar angle distributions. It ranged from 20 to 55 % depending on the 
photon energy. 

In addition to the main background from qq'-j events, a small (less than 5%) fraction 



was due to four-fermion processes which were generated according to EXCALIBUR |T2 
The numbers of selected events and the expected background are listed in Table B 



channel 


(GeV) 


events 


background 


^ ^ 77 


189 


11 


19 ±2 


5 ^ 77 


192 to 202 


19 


24;^ 


S ^ g g 


189 


771 


782 ± 24 


S g g 


192 


113 


113 ±3 


S ^ g g 


196 


339 


316 ±5 


S ^ g g 


200 


342 


330 ±6 


S ^ g g 


202 


169 


158 ±3 



Table 2: Number of selected events for the two decay channels and expected number of 
background events. The background for the 5* — > 77 channel is dominated by the QED 
process e^e~ 77(7), for the S g g channel by the process e'^e^ qq'-f. The errors 
include systematic effects (see text). 



4 Results 

The photon recoil mass spectra obtained for the two decay channels are shown in Fig. |^ 
and Fig. |. The dat superimposed on the expected background distributions. In the 
case of the 5* — > 77 channel, the QED background generator included corrections only to 
order and therefore no additional radiation was simulated. Additional radiation tends 
to give rise to a tail of events having low values of A (Fig. ||). These events were removed 
only from the selected sample of real data, and therefore a corresponding normalisation 
correction of (— 13l7)% "^^^ apphed to the simulated sample. This correction was the 
dominant contribution to the systematic uncertainty for the S* — > 77 channel. 

In the case of the S —>■ gg channel, the systematic error was due to the Monte-Carlo 
statistics and to the uncertainty on the luminosity determination, which was 0.56% for 
the 1998 data and 1.0% for the 1999 data. The e~^e~ qqj background for the 189 GeV 
data was generated with PYTHIA version 5.722, which did not accurately reproduce the 
angular distribution of the radiative photon. Therefore the Monte-Carlo events at that 
energy were corrected on the basis of the ratio between events generated at higher energies 
according to PYTHIA version 5.722 and PYTHIA version 6.125. The systematic error 
on the number of expected events at 189 GeV includes the uncertainty in this correction. 

No excess of events and no clear evidence of anomalous production of events with 
monochromatic photons is observed in either channel. Therefore a limit on the cross 
section of the new physics reaction contributing to the two topologies was set. 

The number of detected events, the background rate and the detection efficiency de- 
pend on the S mass hypothesis considered. In addition, when the expected total width 
for a given ms value is comparable with the experimental resolution or larger, the data 
were compared with the background events in a region corresponding to 80% of the signal 
area. As a consequence, the limit on the signal cross section depends on both ms and 
a/F. To take into account the different sensitivities of the two analysed channels, the 
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Figure 3: Photon recoil mass spectrum for the 777 candidates (points) and the expected 
background (histogram). The average number of entries per event in the data is 2.3. The 
bin size takes into account the experimental mass resolution and the expected signal 
width. 



likelihood ratio method was used ||13[. Since the expected 5* branching ratio and total 
width depend on the mass parameters, as explained above, the 95% Confidence Level 
cross section limit was computed as a function of ms and \fF for the two sets of param- 
eters listed in Table|l|. The result is shown in Fig. ^ By comparing the experimental 
limits with the production cross section computed from (|l]), it is possible to determine a 
95% Confidence Level excluded region of the parameter space. This is shown in Fig. ^. 
As explained in [Q, to keep the particle interpretation the total width F must be much 
smaller than ms and therefore the region with F > 0.5 x 7715 was not considered. The 
95% Confidence Level limits on the cross section times branching ratio for the two decay 
channels are given in Fig. 0. They are obtained for a/F > 500 GeV, corresponding to 
the region where the expected signal width is independent of \/F as it is dominated by 
the experimental resolution. 



5 Conclusions 

The first search for the production of 6*7 (-P7) where S (P) is a CP-even (CP-odd) 
state of the sgoldstino, the goldstino supersymmetric partner, was made using the data 
collected by DELPHI at LEP in 1998 and 1999 at centre-of-mass energies from 189 to 
202 GeV for a total integrated luminosity of about 380 pb~^. The 777 and 'ygg final 
states expected from S{P) 77 and S{P) gg production and decay respectively. 
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Figure 4: Photon recoil mass spectrum for the ''ygg candidates (points) and the expected 
background (histogram). 

were studied. No evidence of a signal was found in either channel. Upper limits on 6*7 
(P7) production in the {ms{mp), a/F) plane were derived. 
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Figure 5: Cross section (pb) upper limit at the 95% Confidence Level as a function of 
ms and \/F for the two sets of parameters of Tab. 0. 
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Figure 6: Exclusion region at the 95% Confidence Level in the ms, \fF plane for the 
two sets of parameters of Tab. |I|. 
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Figure 7: Cross section times branching ratio limits at the 95% Confidence Level for the 
two decay channels investigated. They are obtained for \/F > 500 GeV, corresponding to 
the region where the expected signal width is dominated by the experimental resolution. 
The bin size was chosen to match the experimental mass resolution. 



